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Abstract 

In this Letter we discuss a few issues concerning the magnetic susceptibility of the quark 
condensate and the Son-Yamamoto (SY) anomaly matching equation. It is shown that the 
SY relation in the IR implies a nontrivial interplay between the kinetic and WZW terms in 
the chiral Lagrangian. It is also demonstrated that in a holographic framework an external 
magnetic field triggers mixing between scalar and tensor fields. Accounting for this, one may 
calculate the magnetic susceptibility of the quark condensate to all orders in the magnetic field. 



I. INTRODUCTION 

The magnetic susceptibility x> introduced in [l| in the sum rules framework, is an interesting 
characteristic of the vacuum response to an external magnetic field in the confinement phase. 
It measures the induced tensor current in the QCD vacuum. The expression 
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has been obtained analytically in |2j using the OPE and pion dominance for the (VVA) cor- 
relator of two vector currents and one axial current in the kinematics where two virtualities 
of the external legs are large while one vector current represents the constant external field 
strength. Surprisingly it differs from the sum rule fit by a factor of 3 which implies that some 
qualitative essential effect responsible for the disagreement has been overlooked yet. The sev- 
eral phenomenological estimates yield the low value of the susceptibility while the large N c 
consideration jsj fits the Vainshtein relation ([[]). 

Hence it is natural to look for alternative derivations of x to identify the missed ingredients. 
The problem was discussed in the holographic hard wall model involving the 5D Yang-Mills 
and Chern-Simons (CS) terms. The Lagrangian corresponds to the gauge theory on the flavor 
branes extended along the radial coordinate in the AdS space. It turns out that in this model 
the Vainshtein relation ([T|) is not exact, however it is fulfilled with good accuracy j^]. Moreover 
it was shown in 4j that the whole answer follows from the CS term which implies that we are 
dealing with a sort of "anomalous" phenomena. 

Other more refined holographic models have been considered by Son and Yamamoto (SY) 
in j^J. They derived a new relation between two-point and three-point correlators which yields 
nontrivial matching conditions for the low-energy QCD parameters of the mesons. The SY 
relation is assumed to be valid at any momentum transfer, for instance, Vainstein expression 
(PQ) follows from the SY relation at large virtualities if one assumes that operator product 
expansion of QCD is applicable. This 'if is important because the SY model does not support 
OPE per se: dependence on momentum transfer is exponential and does not contain power 
terms required by OPE. Different aspects of the SY relation were discussed in Q-8]. 

The situation looks a little bit puzzling since there is no field theory derivation of the SY 
relation yet. The expression for \ m terms of pion decay coupling suggests that it can be 



obtained purely in terms of the chiral Lagrangian together with the SY relation. With the 
holographic experience it could be expected that the Wess-Zumino-Witten (WZW) term in the 
chiral Lagrangian related to the CS term in 5D should be responsible for the nontrivial answer. 



In this note motivated by the comments above we consider the additional arguments con- 
cerning the derivation of x- Since the holographic model of QCD is nothing but the extended 
chiral Lagrangian it is natural to look more carefully at the place of the SY relation in the 
ChPT per se. The small Q 2 region is the most comfortable to be analyzed in ChPT hence we 
shall look at the first terms in small Q 2 expansion. The relation between the ChPT parameters 
at the tree level has been discussed in jf]] and we extend it to the one loop level focusing at the 
chiral logs. It will be shown that the SY relation holds true for the log terms. 



A simple argument involving the calculation of the quark determinant in the tensor source 
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background implies that the nonvanishing magnetic susceptibility corresponds to the pecu 
additional mixed term in the Lagrangian. In the improved holographic model for QCD {9- 
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the tensor source in D = 4 theory is promoted into the tensor field in D = 5. We shall analyze 
the improved model in a magnetic field focusing at the scalar-tensor mixing. It turned out that 
in the improved model the magnetic susceptibility can be obtained to all orders in the magnetic 
field. In a small field the magnetization grows linearly with the field, in accordance with its 
generally established properties, while in a large field it does not depend on the magnetic field. 



The paper is organized as follows. In Section 2 we consider the SY relation in ChPT and 
show the matching of the chiral logs in this relation. Some general comments concerning the 
SY relation are also presented. In Section 3 we consider the improved holographic model for 
QCD with the tensor field and demonstrate how the scalar-tensor mixing in the magnetic field 
yields the magnetic susceptibility. Some discussion concerning the proper degrees of freedom 
can be found in the last Section. 
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II. SY RELATION WITHIN THE CHIRAL LAGRANGIAN 



A. SY relation and chiral logs 

Let us analyze the SY relations in the framework of ChPT. The SY relation has been 
obtained holographically and is based on the simple D = 5 action on the worldvolume of probe 
flavor branes involving Yang-Mills and Chern-Simons terms for the flavor gauge group. Taking 
into account that correlators of the vector and axial currents are two independent solutions to 
the second order differential operator in the radial coordinate z in AdS space their Wronskian 
is ^-independent. This argument works when the CS term is neglected. On the other hand CS 
term itself yields the nontrivial (VA) correlator in the magnetic field which is proportional to 
the same Wronskian. Hence in a weak magnetic field the following relation holds [5( 

™ T (Q 2 ) = ^ - ^ [n A (Q 2 ) - u v (q 2 )] (2) 

where wt is defined via the two-point correlator in the external weak electromagnetic field with 
the constant field strength F a/3 

(VpAv) F = ~^ M? 2 )(-<Z 2 ^ + QvtfFp* ~ %<fh*) + w L (q 2 ) q v q°F^\ , (3) 

V,A are the vector, qVy^q, and axial, q Ay^q, currents, F denotes the dual field strength, 
F-ys = \ t-fSapF - 13 ', and Ua , Hy are the corresponding two-point correlators. The relation holds 
for all values of Q 2 = —q 2 . 

The definitions of the above correlators are 

l -Tx{QVA){V,Au) F = j d A xe t<lx (T{V^x)A u (0)}) P , 

l -Ti(yV)Ii v {Q 2 ){ qil q u -g^q 2 ) = j d A x e^ x (T{V,(x)V u (0)})o , (4) 

l - Tr (AA) U A (Q 2 )( qtl q u - 9llu q 2 ) = J d A x e igx (T{A,(x)A u (0)}) , 

where flavor dependence on matrices of vector and axial currents, V and A, as well as that for 
the electric charge, Q, is factored out. 

It is natural to look at the matching of the SY relation with the ChPT since the holography 
provides the ChPT derivation from the "first principles" . The SY relation is derived from the 



D — 5 equations of motion hence according the holographic dictionary it should correspond 
to the tree approximation in the ChPT. The condition imposed by the SY at Q 2 = on the 
parameters of the chiral lagrangian has been found in |f| 

L 10 = -4tt 2 C 22 (5) 

where L w corresponds to the even term in the chiral Lagrangian at the 0(p 4 ) order while 
C22 corresponds to a particular odd term at the 0(p 6 ) order. The condition is unexpected 
since it relates the odd and even terms in the Lagrangian. Unfortunately this relation between 
constants can not be used as the test of the SY relation since C22 is not known with the high 
accuracy. 

To get some test of the SY relation let us focus at the terms in the correlators involving 
the chiral logs. Some comments are required before the looking at log terms. Naively such 
terms are subleading in 1/N C hence in the holographic approach these should be considered as 
corrections to the equations of motion. On the other hand the log terms are considered as the 
renormalization of the constants in the chiral Lagrangian hence one should assume that the 
relation (jSJ) valid at the tree level holds upon the renormalization. Therefore it is natural to 
look at the matching of the chiral logs. 

At the right hand side of Eq. ([2J the chiral log follows from the pion loop in the correlator 
of the vector currents 

< ir (Q 2 ^0)=clogQ 2 , c=-^. (6) 

There are no logs in the correlator of the axial currents. On the other hand the chiral log in 
the (VVA) correlator can be traced from the particular term in the WZW term in the chiral 
Lagrangian responsible for the decay 7* — > 3n, 

S'wzw = / TrA(dU-'U) 3 -J^ J d 4 x Tr F^r «V d 5 n . (7) 

Converting two pions from this vertex to the vector current and associating the remaining pion 
with the axial current we get the log Q 2 contribution to wt , 

< ir = Cl logQ 2 , Ci = ^c, (8) 
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which is consistent with the SY relation. 

Let us emphasize that there is no freedom in the terms in the chiral lagrangian involved into 
the chiral logs, hence the matching is exact although at the subleading order in 1/N C . Note 
that similarly to the tree-level case, the SY relation implies an unexpected relation between 
the coefficient in front of the even kinetic term at the 0(p 2 ) order and of an odd WZW term 
at the 0(p A ) order in the chiral Lagrangian. 

B. The mixed term in the chiral Lagrangian 

Let us argue that nonvanishing magnetic susceptibility implies a peculiar term in the effective 
Lagrangian. To this aim we introduce a source term for the quark tensor current into the QCD 
Lagrangian 

8Lt = B^qBa^q = iB^qBa^q (9) 

where B^ v is an external source field whose possible interpretation shall be discussed below 
and B is a diagonal flavor matrix. Accounting for the chiral features of the quark operator in 
Eq. ([9]) it simple to determine the corresponding term in the chiral Lagrangian in the linear 
approximation in the B^ field, 

SL WZW = ~X (ft) B, V F^ Tr (U + U*)BQ , (10) 

where x is the magnetic susceptibility, (qq) is the quark condensate and U = exp(2in a t a / f n ) 
is the mesonic matrix (f w = 92MeV). This can be viewed as a definition of the magnetic 
susceptibility. Note that this term to some extent can be considered as the shift of the effective 
quark mass in the external fields. In the next Section we shall see that this term promoted into 
the holographic D = 5 action provides an important scalar-tensor mixing. 

It it worth making a few comments concerning the implications of this effective WZW-like 
term in the chiral Lagrangian. First, the vacuum tensor current proportional to the chiral 
condensate in the magnetic field can be attributed to the stringy degrees of freedom if we 
identify the tensor source in f lTUj) as NS or Ramond two-form fields. With such identification 
the flow of the Fl or Dl strings in the vacuum occurs in the magnetic field. On the other hands 
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the mesons are identified hole-graphically as the Fl strings connecting the flavor branes hence 
such "stringy" current corresponds in fact to a kind of mesonic vacuum current. 

Secondly, there is an anomalous electromagnetic current proportional to the condensate in 
the external tensor field. Indeed we defined the current 

= ^ (id 

which gets contributions from the sources of the tensor field or from the varying pion field. 

Wb = \x {m) & [B^Tr (U + rf)BQ] . (12) 

For the varying pion case the anomalous current is the analogue of the Goldstone-Wilczek cur- 
rent. For the varying tensor field an interesting possibility emerges. Using the relation between 
the massive vector and tensor in four dimensions we could get the nonvanishing electromagnetic 
current if the vector meson gets condensed. There are some indications of such a condensation 



in the magnetic field both in the effective theory 12j] and in the holographic framework 13 1. 
Hence one could speculate about the nonperturbative current proportional to the product of 
quark and vector meson condensates. 

Finally, if we expand the anomalous term in the pion field we could get the anomalous 
interaction of pions with the tensor current in the magnetic field. For instance, the matrix 
element 

(0|ga^g|7r 7r°) = X (qq)F^ . (13) 
C. On the derivation of the SY relation 

The SY relation has been obtained in a slightly tricky way, hence it would be nice to get it 
more regularly as a kind of a Ward identity. Here we restrict ourselves by two generic remarks. 
Since the key observation in the derivation in j^] was the z-invariance of the Wronskian of the 
vector and axial currents it is reasonable to look at the radial variable z in the Hamiltonian 
framework. That is, following [l^| we assume that it is considered as a time variable for the 
RG Hamiltonian evolution in the D = 5 gauge theory. 

In the Hamiltonian framework of the gauge theories there are two natural equations involv- 
ing the dependence on the boundary values of the dynamical variables. These are the gauge 

7 



constraint or Gauss law and the Hamiltonian constraint or a kind of the Hamilton- Jacobi (HJ) 
equation. We are in a peculiar situation with the Hamiltonian constraint since the metric 
depends on the radial coordinate and is therefore "time-dependent" . 

First consider the Gauss law constraint with respect to the flavor gauge group SUl (Np) x 
SU/?(iVV) on the flavor branes. In the Hamiltonian approach the Gauss law reflect the gauge 
invariance with respect to the flavor gauge group and can be identified with the generator of 
the z-independent gauge transformations. Since there are D = 5 CS terms for the left and right 
gauge flavor fields the canonical momenta get modified 



while the canonical momenta for the scalar field are standard. Using the Hamiltonian relation 

n = £, (15) 

OX 

where S is the action, and the standard holographic relation for the D = 4 currents = 44 
one immediately recognizes that the Gauss law constraint in the bulk theory precisely produces 
the anomaly equation for the axial current at the boundary including the mass term. The fact 
that the Gauss law is valid at any time in D = 5 theory gets translated into the claim that the 
axial anomaly is seen at all scales in the boundary D = 4 theory. Note that in the conventional 
gauge theory the Gauss law is complemented by the gauge A = 0. In the current situation 
the similar equation reads as A z = 0; however, one should not forget that the pion field can be 
identified with the holonomy of the radial component of the flavor gauge group. 

In the holographic setting the HJ equation for the bulk metric has been identified as the 

n 

RG equation in the boundary theory in [14J]. Here we have to consider a similar HJ equation 
for the gauge fields and scalars. Taking into account the shift of the canonical momenta and 
forgetting for a moment the metric one obtains for the left gauge part of the total Hamiltonian 



and similarly for the contribution of the right gauge field and scalars. The HJ-type equations 
are quite convenient for the derivation of Ward identities in the boundary theory since it 
involves the desired variational derivaties. It is important that the HJ-like equations due to the 



Ii L = E L + A L F L , 



T1 R = E R - A L F L , 



(14) 




(16) 
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change of the canonical momenta involve the terms with the different number of the variational 
derivatives. Hence potentially one could hope that the additional variational derivatives of the 
HJ equation upon taking into account the Gauss law constraint would yield the SY relation. We 
did not succeed along this way of reasoning, however we plan to discuss the complete set of the 
Ward identities induced from the bulk theory elsewhere. In particular we plan to elaborate the 
constraints emerged from the dynamics of the higher rank fields induced by the color branes in 
the brane approach. Note that some examples of the derivation of the boundary Ward identities 
from the bulk HJ equation can be found in [jjj. 

There is also some analogy with the N = 1/2 SUSY YM case which can be considered as 
N — 1 SYM theory in the self-dual constant graviphoton background C\. The following term 
gets induced in the graviphoton field 

5L = dd A FA A (17) 

which is analogous to the anomalous term in QCD we have discussed. The analogy with QCD 
becomes even more close when we remind that the gluino condensate is developed in N = 1/2 
SUSY YM like the chiral condensate in QCD. Moreover in the N = 1/2 theory one can consider 



This Ward identity amounts to a 



17] . Since the spectral 



the Ward identities reflecting the single unbroken SUSY 
particular degeneration in the spectral densities in the J = 1 ± channels 
densities follow from the two-point correlators these Ward identities can be considered as some 
analogue of the SY relation in QCD without the anomalous three-point correlator. 

III. A HOLOGRAPHIC MODEL WITH THE TENSOR FIELD 

In this Section we shall consider the scalar-tensor mixing in the improved holographic model 



of QCD which involves the tensor field 



-|ll|. It is a 5-dimensional gauge theory embedded in 



a pure AdS geometry with an infrared hard-wall boundary: 

I 2 

ds 2 = — (-dz 2 + rj^dx^dx") , < z < z m , (18) 

where rj^u is mostly negative: t] = diag(H ), and I is the AdS^ radius and shall be omitted 

hereforward (thus rescaling the coupling constants). This model contains three types of fields: 
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a complex scalar X, two gauge fields L M and R^, and a complex antisymmetric tensor B^ v . 
They are put into correspondence with the following operators of QCD: 

Qr } B lu I ' 9£ <?7/^{ <"> L i s > ( 19 ) 

where a^ u = |[7 M , 7^]; /, / are the flavor indices of QCD with respect to the (global) U(Nf) l x 
U(Nf)n symmetry which becomes the gauge group of the five-dimensional theory. Accordingly, 
the fields X and B^ u are bifundamentals, whereas and are adjoint with respect to 
U(Nf) L and U(Nf) R . These properties allow us to properly introduce covariant derivatives: 
DX = dX - iLX + iXR, H = DB = dB - iL A B + iB A R. 



The action proposed in llj is: 

- 2g B Q (BAH + -B + AH)+ m B \B\ 2 ^ + ^ (X + F L B + + c.c.) j . (20) 



If 



which takes into account the 



This action is a modification of a simpler hard-wall action 
tensor field. The interaction term XFB on the AdS boundary is reduced to the term (]10p of 
the chiral Lagrangian which we have discussed in the previous Section. 

The constants have been fixed in previous works by comparing various correlators at large 



Euclidean Q 2 with OPE in QCD [ll|, ll8H20j . The masses are fixed by requiring that the scaling 
properties of the fields match those of the corresponding operators in the UV: m 2 x = —3, 
rriB = 2. Note that due to a non-canonical form of the kinetic term of the tensor field its 
physical mass is actually 1 in units of i~ x . In this case the vacuum solution for X is 

X(z) = \ (mz + 4- (qq) z 3 ) x l NfXNf . (21) 
^ v 9x 

From now on we shall only consider the Abelian degrees of freedom, as the flavor structure of 
the 5D fields is trivial, since the condensates, both scalar and tensor, as well as the electroma- 
gentic field, are diagonal in the flavor space. Hence, the equations of motion for each individual 
flavor qf are the same as for the singlet component with a substitution Fff^ — > ejF^ 1 ^ , where 
Cf is the electromagnetic charge of g/. Furthermore, we shall be working in an ansatz where 
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the axial field [Lm — Rm)/V^ is zero, which is consistent with the equations of motion. We 
are left with only the vector field Vm = (L M + R M )/y/2. In this case the covariant derivatives 



become ordinary. Let us also split the scalar and tensor fields into real and imaginary parts: 
in QCD: 



X = — + ~^ ; Bmn — - — + ""^/^ — ' These new fields are dual to the following operators 



1 

i 

q^q^V^. (22) 



qq<r>X+, -j= qa^q «-> B +fXV , 

i 

iqi$q O X_ , — gjsa^q B_^ u , 



We have noted that is bifundamental with respect to U(Nf)L x U(Nf)n which guaran- 
tees its being complex-valued, its real and imaginary parts corresponding to the tensor and 
pseudotensor operators (|22|) . These operators happen to be related to each other in 4D, 
qv^l^q — -^ V \ p q aXp q^ that fact is reflected in Eq. ([9]). From the holographic point of view, 
this condition is ensured by the fact that the kinetic term for (|20|) is of the first order in 



derivatives, which leads to its complex self-duality 10 



llj . Thus, as we shall see, the "double 



counting" of the degrees of freedom that arises after we have introduced a complex tensor field 
is compensated by constraints imposed on half of them, see Eqs. ( 12 7p . One may wonder whether 
it is possible to avoid this redundancy by dealing with a real-valued tensor field from the be- 
ginning; however, amending the model in this way while preserving holographic field-operator 
correspondence rules and general self-consistency appears to be quite cumbersome. 

Let us now rewrite the action ( 12"U|) in terms of those fields: 



S = J d 5 x^Tr{-^F* + 9 fe MNp Q R (B„ MN H, 



PQR — B + mnH-PQr) 

A 



[~ 9 B m B B ±MN B ±MN + 9 -f {d M X ± d M X ± - m\Xl) + ±X ± [F V ) MN B^ N ] }. (23) 
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A. Equations of motion 



The action (1231) yields the following first-order equations of motion for the tensor field: 

± ze^ Xp H ±zXp + 25f = -^X^Ff, 

± z e ^ P ° H + 2 Bg = ———XzfFy Z , (24) 
3 + Ag B 

where the indices are contracted with a flat metric diag(H ). They may be rewritten 

as second-order equations in which the real and imaginary components are disentangled and 
(Fy)^ is assumed to be z-independent (which is a self-consistent solution in the case of a 
constant uniform magnetic field). Along with the equation on X we have: 

zd z (zH z ± a ?) + Bf + z%H^ = -A- \x ± Ff ± zd z X T F v 



A 



z 2 d x H^ z + B£ = ±^-zd x (X T F» X , , 
*9b 



d z ( \d z x±) + A x± - -a,d»x ± = -WfvUb^. 

\z 6 J z b z 6 g x z 



(25) 



Directing the third axis along the magnetic field so that (F v ) 12 = (F v ) 03 = B we get the 
following equations on (B±) 12 and X±: 

(d 2 z + -d M -\- d,dA (B ± ) l2 = ~\x ± (F v ) 12 , 
\ z z l J ogs z 

(dl --d z + - 2 - d,dA X ± = (F v ) 12 (B ± ) 12 . (26) 

\ z z j g x 

From Eqs. ( 125]) it follows that we will also have nontrivial (B T ) 03 , (B T ) 0z , and (B T ) 3z compo- 
nents, which may be expressed through (B±) l2 with the use of Eqs. ( 124"]) (assuming e 0123 = 1): 

(£ T ) 03 = ±zd z (B ± ) 12 , (B T ) 0z = ±zd 3 (B ± ) 12 , (B T ). 3z = ±zd (B ± ) 12 . (27) 

A most general property of the equations is that the scalar and tensor degrees of freedom 
X + ,B +12 decouple from the pseudoscalar and pseudotensor X_,£L 12 , thus forming two inde- 
pendent sectors, while due to complex self-duality B t03 , B T oz, B t3z are admixed to the first 
(second) sector of the solution. Those sectors may be treated independently. 
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B. Solutions and boundary conditions 



After we Fourier-transform the equations, solutions of the Eqs. f !2B| l2T)j) assume the form: 

X + + iX_ = z 2 f x {qz)e ikx *- iu)t - 
(B + + iB_) u = -^f B ( qz )e^-^, (21 



where fx(qz) an d fs(q z ) are, generally speaking, superpositions of four Bessel functions: 



f x (qz) = dJ^y/T+Pqz) + C 2 J x (s/l^qz) + C^y/l + fiqz) + C^y/T^qz); 
f B (qz) = CiJ^y/T+Pqz) - C 2 Ji(y/l - Pqz) + C 3 yi(y/T+Pqz) - C^y/T^qz) , (29) 
where (3 = ■ — ■ |B/g 2 | and q 2 = u 2 — k 2 is the Minkowski 4-momentum squared. J , y 



are the Bessel and Neumann functions or their analytical continuations (if we want to consider 
greater magnetic fields or solutions with Euclidean momenta). The Neumann functions in 
( l2"9l correspond to non-normalizable modes of the AdS§ fields ( 1281) in the UV. As we are only 
interested in the mixing between vacuum expectation values without sources, C3 = C4 = 0. C\ 
and C2 are determined by the boundary conditions in the IR. 

According to f l27|) . there is a constraint that relates (5 =F ) 03 to (B±) 12 , which means that in 
order to construct a self-consistent variation principle for the tensor field one needs to take into 
account that half of the tensor degrees of freedom are not independent due to the tensor field's 



complex self-duality. Such variation principle has been proposed in ll|, and it states that in 
our case 

5 B S = 2g B J d 4 x (B +12 + £L 03 ) 5 (B 12 - £L 03 ) . (30) 

Note that according to f )27|) -B+12 an d -B-03 have equal normalizable modes, thus contributing 
equally to the tensor condensate (1301 . Since the kinetic term of the tensor field is of the first 
order in derivatives, the boundary variation term in ( |30l) contains no differentiation with respect 
to z. Hence it makes sense to impose on it a Dirichlet boundary condition at z = z m rather 
than a Neumann one. From ( 1301) it also follows that a Dirichlet condition has to be imposed 
on the sum (-8+12 + -8-03)- Hence, 

Ci = MVT=P qz m ) + y/T^g qz m Jl(y/T=V qz m ) 
C 2 Ji(vTT? qz m ) + ^TT? qz m Jl\y/T+Vqz m ) ' 
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There is no infrared boundary condition for X, so the overall value of Cj remains undeter- 
mined. Nevertheless, we can obtain the ratio of the tensor and scalar condensates (the scalar 
one is determined from (}2~Tj) . while the tensor condensate is read off of the variation of the 
action with respect to the tensor field (130]) ): 



9x f B (qz) 



(qa 12 q) oc 8g B - lim 



(??) oc g\ lim 



fx{qz) 



z->0 z 



hence 



MB;g) = i^ = ^lim /s( ^ B) 



(qq) 9x fx(qz; B) ' 



(32) 
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FIG. 1. Magnetization /u(B) (blue) vs its strong field asymptotics (red). 



Setting the 4-momentum to zero we are able to obtain the magnetization /i(B) and the 
magnetic susceptibility x(B) = -=r yLt(B) for a uniform condensate in terms of Bessel functions 



of 



-Hz^- They are presented here on Figs.[T]and[2l respectively. 
Their main properties are that 



lim /u(B) = — V3b 

B^oo g x 



(33) 
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FIG. 2. Magnetic susceptibility of the quark condensate x(B) 



and 



X(B) 



1 A 2 



-b 2 4 + o (b 4 4)) , 

• J 



B -> 0. 



(34) 



ff|4 V 96 QxQb 

Note that the magnetization changes its behavior from one linear in the magnetic field to a 
constant at values of the magnetic field of order of B ~ z~ 2 ~ Aq CD . Its constant asymptotic 
is a behavior to be expected. At large magnetic fields the dynamics of the theory become 
effectively two-dimensional and the tensor chiral condensate is kinematically reduced to a scalar 
one. 

If we substitute the values of the constants of the model from [ll| we obtaii|j] 



At(oo) = 1/V3 . 



X(0) 



m 

72 



Let us recall that z m is fixed by the mass of the p-meson 



that x(0) 



-0.08 m~ 2 ~ 



0.13 GeV . One may compare that to the results 



(35) 

2.4 m" 1 , which means 



2l|, where the 



susceptibility has been analyzed from the point of view of sum rules and has been determined 



1 Our constants gs and A are 3 times larger than those in 



ll| due to the fact that our variation of the action 



pop with respect to the tensor field is, in a similar way, 3 times larger than the variation in 



15 



to be y ~ —3.15 ± 0.30 GeV -2 , while the pion dominance and OPE for the (VVA) diagram 
give [2J X ~ — 8.9 GeV" 2 . (Other results include a holographic computation of the (VVA) 
diagram, which led to a value x ~ —11.5 GeV -2 J4|; the use of vector dominance 22 j gives 
X ~ —(3.38^-5.67) GeV~ 2 .) One can see that there is a large discrepancy between the numerical 
results that clearly requires more investigation. However and more importantly, our results 
reproduce the general properties both of the susceptibility and of the magnetization - the 
weak-field expansion of the former and the negative constant asymptotic of the latter. 



As for the particular value o 
saturation value in a large field 



the magnetization, the lattice calculation yields a different 
231 ] : lim //(B) = — 1. It has also been discussed in the NJL 

B— >oo 

model 24J. In theory, such a value would tell us that both condensates get contributions only 
from the LLL (lowest Landau level) in the strong magnetic field. On the other hand a different 
saturation value obtained in this paper implies that the picture is more complicated. Indeed, 
as it has been shown by Miransky et al. 25[, in some problems the summation of the infinite 
number of hig her Landau levels is needed to reproduce the correct result. It has also been 
argued in 



that the LLL approximation is reliable only in the kinematic region when the 
momenta satisfy the condition q\ 3> q». In our particular case the problem is completely 
static (qi = 0), so this condition is not fulfilled and the careful treatment of the higher Landau 



levels is desirable. Another simp 
the Dirac operator spectrum 
the result obtained in 



e argument concerns the derivation of the magnetization via 



231 ] . An analogue of the Casher-Banks formula implies that 



23] is based on the factorization of the product of two operators under 



the averaging over the gluon configuration. The lack of factorization could be the origin of 
disagreement with the lattice result. Anyway, this point needs further clarification. 
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C. Vector current 



Let us now consider perturbations of the vector field about the solution Fyn = B. They 
obey the following equations: 



d, I -d z V 3 



-d 2 V 3 + -d Q d 3 V = 4g 2 \^ 

z z 



( -x 

z 



±-"±32 



-d z 



-d z Vr 



-d 2 3 V 



: d d 3 V 3 = 4g 2 \J2 



d z -X±B± 02 



-do (X±B± 03 ) 



-do (X±B± 03 ) 



• (36) 



One may note that the longitudinal components B± z/1 (as well as B±q 3 ) become sources for the 
vector current and charge density However, the r.h.s. of the Eqs. (I36|36p contains products 



of fields from different sectors of the solution. Furthermore, if we consider small fluctuations, 
the vector field turns out to be a fluctuation of the second order. There are obvious similarities 
with Eq. (Tl2|) . where ir = arg (X + + iXS). 



IV. DISCUSSION 



In this paper we have discussed a few issues concerning the magnetic susceptibility of the 
quark condensate. We have shown that the SY relation which yields the value of the suscepti- 
bility at large Q 2 is consistent with the chiral log counting at small Q 2 . The nonvanishing value 
of the susceptibility implies a specific term in the effective lagrangian and we have analyzed the 
role of this term in the holographic approach. It turns out that it provides the magnetization 
at any value of the magnetic field. Surprisingly the saturation value in the strong magnetic 
field is small and disagrees with the lattice simulation. This disagreement deserves additional 
studies. 

A satisfactory explanation of the SY relation is still absent. It implies a peculiar relation 
between the kinetic and topological terms. Such relation is natural from the brane viewpoint 
and one could expect a kind of Ward identity to stand behind it. We have not found the 
symmetry which would provide such a Ward identity, however more efforts could be made in 
this direction and we plan to return to this point elsewhere. 

Which vacuum excitations are responsible for the magnetic susceptibility? This question 
can be rephrased as one concerning the localization of the quarks involved into the composite 
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operator on some vacuum defects excited by the external magnetic field. The answer potentially 
depends on the interpretation of the background C^ v field. There could be several interpreta- 
tions. If it is the two-form field in NS or RR sectors it would mean that the Fl or Dl degrees 
of freedom are under the carpet. The variant with NS B^ v field has some trouble since in this 
case we are dealing with noncommutative field theory and the field enters other terms in the 
Lagrangian. Hence it is unclear if it would be possible to separate the desired term in the 
effective action in a clear-cut way. If we choose the RR two-form field C2 the product C2 A F 
follows from the CS term immediately, however the proportionality to the chiral condensate 
needs an explanation. 

Finally let us comment on the possible interpretation of the tensor source as the curvature 
of the graviphoton one-form RR field C\. C = dC\. The degrees of freedom naturally charged 
with respect to the graviphoton field are the DO branes, hence one could question how DO 
particles or instantons are captured by the magnetic field. The potential object which could 
be relevant is the dyonic instanton, that is, a blown up instanton with a string attached to it. 
Upon the blow up it behaves as a magnetic dipole with a topological charge. 
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